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—  V*  Combustion  efficiency  of  aluminized  propellants  in  solid  rocket 

motors  is  reduced  by  incomplete  aluminum  combustion  and  two-phase  nozzle 

flow  losses.  'tIic  ^orabusclon  of  klurainize'd*  propellants  can  produce  large 

-  AI/AI2O2 ^agglomerates.,  The  agglomerate  dynamics  within  the  combustion 

chamber  and  nozzle  have  a  significant  influence  on  the  overall  corabus- 

i  —5 

tlon  efficiency  of  the  motor_  .  Agglomerates  are  subjected  to  large 
aerodynamic  loads  within  the  rocket  nozzle  where  the  gas  phase  experi¬ 
ences  a  more  rapid  acceleration  chan  the  agglomerates.  The  drag  load 
deforms  the  agglomerates  and  If  of  sufficient  magnitude  result  in 
breakup.  The  smaller  fragments  have  faster  velocity  and  thermal  equi¬ 
librium  tiroes  and  have  higher  combustion  rates.  For  maximum  combustion 
efficiency  the  aluminum  fragments  must  completely  oxldlzei  and  achieve 
kinetic  and  thermal  equilibrium  with  the  gas  phase.  As  a  direct  result 
of  agglomerate  breakup,  Che  aluminum  combustion  rate  is  increased,  and 
the  thermal  energy  released  is  more  efficiently  transferred  into  exhaust 
kinetic  energy.  -  >  '  :  — 

Photographic  observations  obtained  in  windowed  rocket  raotors^”^ 

8—  1 0 

and  combustion  bombs  indicate  relatively  large  agglomerates  (lOOpm  - 
500pm)  are  formed  on  the  propellant  surface  and  entrained  in  the  com¬ 
bustion  flow;  however,  particle  size  measurements  obtained  from  sampling 
the  exhaust  plumes^ ^  indicate  small  mean  particle  diameters 
(<10uiu).  These  small  exhaust  plume  particles  apparently  result  from  the 
breakup  of  the  larger  agglomerates  during  the  nozzle  expansion 
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process.  Observations  of  agglomerate  breakup  in  a  laboratory  scale 
rocket  nozzle  revealed  an  adequate  correlation  with  Weber  number; 
however,  neither  the  physical  process  of  breakup  nor  the  fragment  size 
distribution  was  resolved^^. 

The  objective  of^this  research  its  to  obtain  physical  data  to 
characterize  the  mechanisms  of  aerodynamic  droplet  breakup.  Experiments 
have  been  completed  In  which  conventional  liquids , (water ,  alcohol  and 
glycerine/water  mixtures)  "and  a  liquid  metal  (mercury)  was  studied.  The 
primary  goal  of  the  conventional  liquid  experiments  was  to  examine  the 
effect  of  liquid  properties  (viscosity  and  surface  tension)  on  the 
breakup  mechanism,  time  scale,  and  fragment  size  distribution.  The  goal 
of  Che  liquid  metal/  experiments  was  to  examine  the  effect  of  the  much 
higher  surface  tension  more  characteristic  of  liquid  aluminum. 

A  key  element  of  the  experimental  effort  is  the  use  of  - 
nonintruslve  laser  diagnostics  including  pulsed  laser  holography ; and 

ni.v'j. 

laser  Doppler  veiociraetry^^  The  exceptional  temporal  and  spatial  resolu¬ 
tion  of  pulsed  laser  holograph^  provided  the  ability  to  resolve  the 
mechanism  of  br^kup  and  the  size  distribution  of  the  fragments,  ^aser 
-Uoppler  velocimetrjr  was  used  to  determine  drop  velocity  distributions 
along  the  nozzle  revealing  the  rapid  acceleration  of  the  flattened 
droplets  and  then^ surprisingly^  the  milder  acceleration  of  the 
fragments.  ^  The  high  drag  to  mass  ratio  of  the  flattened  droplets  was 

I 

expected  aijd  resolved  during  the  initial  phase  of  the  breakup  process. 

I 
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2.0 


EXPERIMENT  CONSIDERATIONS 


Aerodynamic  droplet  breakup  is  characterized  by  the  magnitude  and 
duration  of  the  aerodynamic  forces.  For  nozzle  contractions,  the 
droplets  are  loaded  as  a  result  of  the  higher  gas  phase  acceleration. 

The  droplets  respond  to  the  aerodynamic  forces  by  deforming  and 
accelerating.  The  initial  acceleration  scales  with  the  initial 
diameter:  however,  as  the  pressure  forces  flatten  the  droplet, 
increasing  its  cross-sectional  area,  the  drag  load  and  the  acceleration 
increases  rapidly.  This  sequence  of  events  can,  of  course,  lead  to 
catastrophic  deformation  and  droplet  breakup  if  the  slip  velocity  is  of 
sufficient  magnitude  and  duration. 

The  liquid  surface  tension  is  used  to  scale  the  aerodynamic  load 
forming  the  Weber  number: 


We 


p  Av 


(1) 


p  =  density 
Av  “  slip  velocity 
(|)  =  diameter 
o  “  surface  tension 
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The  duration  of  the  aerodynamic  loading  can  be  scaled  by  the  natural 
period  for  hydrodynamic  oscillation,  which  is  a  measure  of  the  droplet 
response  time: 


T 

n 


Pd  ♦d 


3  1  jZ 


(2) 


For  nozzle  accelerations,  an  additional  time  scale  exists;  namely  the 
time  required  for  velocity  equilibrium,  Assuming  Che  droplet 

Reynolds  number  relative  to  the  gas  phase  is  sufficiently  high  (l.e., 
Re  >  10^  ),  the  drag  coefficient  for  rigid  droplets  is  approximately 
unity  and  Che  corresponding  time  scale  is 

"d  ♦d  ( 

V  pAv 


However,  since  the  liquid  droplets  are  not  rigid,  the  droplet  cross 
section  and  the  associated  drag  load  Increase  rapidly  within  the  nozzle 
contraction.  Therefore,  the  time  scale  for  liquid  droplet  velocity 
equilibrium  is  no  doubt  faster  than  the  rigid  droplet  estimate. 

A  few  liquid  types  were  selected  to  examine  the  effect  of  liquid 
properties  such  as  surface  tension  and  viscosity.  For  the  conventional 
liquids,  water  and  alcohol  were  selected  for  the  surface  tension  series; 
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and  water  and  a  glycerine  water  (57/43)  mixture  were  selected  for  the 
viscosity  series.  The  mixture  was  selected  to  achieve  a  viscosity  ten 
times  that  of  water. 

A  second  series  of  experiments  were  designed  to  examine  the 
breakup  dynamics  of  high  surface  tension  liquid  metals.  Mercury  was 
selected  because  of  its  compatibility  with  the  droplet  generation 
technique  and  wltn  room  temperature  air  flow. 

2.1  Nozzle  Facility 

An  e;'periraent  was  designed  in  which  droplet  breakup  could  be 
investigated  in  a  nozzle  contraction.  The  nozzle  was  designed  to  allow 
droplet  Injection  upstream  of  the  contraction.  Windows  formed  the  sides 
of  the  two-dimensional  nozzle  allowing  observation  of  the  droplet  inter¬ 
action.  The  gas  was  delivered  from  a  settling  chamber  through  a 
sintered  metal  throttle  to  the  nozzle  plenum  where  the  droplets  were 
injected.  The  nozzle  area  ratio  was  eight.  Figure  1.  For  the  conven¬ 
tional  liquid  experiments,  the  nozzle  was  operated  at  subsonic  exit 
velocities.  However,  for  the  liquid  metal  experiments,  higher  dynamic 
pressures  were  required  to  achieve  droplet  breakup.  Higher  droplet 
loading  was  achieved  by  increasing  the  gas  velocity  and  hence  the  slip 
velocity.  In  fact,  the  nozzle  exhaust  was  modified  to  allow  sonic  flow 
in  the  nozzle  throat  for  the  liquid  metal  experiments,  although  droplet 
breakup  occurred  upstream  of  the  throat. 

Piezoelectric  droplet  generators  were  constructed  to  produce 
highly  monodlsperse  droplet  streams.  The  basic  design  was  that  of 
modern  ink-jet  printers.  The  complete  droplet  generator  is  composed  of 
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Figure  1.  Nozzle  Facility  with  Piezoelectric  Droplet  Injector 


a  liquid  supply  tank,  an  electrical  driver,  and  the  ink-jet  devices 
(Figure  1).  The  size  of  the  droplets  generated  by  our  devices  was 
variable  from  about  1.5  to  3.0  times  the  orifice  diameter  by  decreasing 
the  excitation  frequency.  For  conventional  liquids,  droplets  were 
generated  from  160  pm  to  400  pm  in  diameter,  and  for  mercury,  droplets 
were  generated  from  190  pm  to  560  pm.  The  droplet  generator  operating 
conditions  are  listed  in  Table  I  for  the  conventional  liquids  and 
mercury. 


TABLt;  I.  DROPLET  GENERATOR  OPERATING  CONDITIONS 


DIAMETER 

EXCITATION 

LIQUID 

TIP 

(urn) 

DROPLET 

(pm) 

FREQUENCY 

(kHz) 

Water 

100 

160 

10.6 

120 

220 

8.0 

190 

410 

0.8 

Alcohol 

100 

164 

11.2 

150 

265 

9.0 

Glycerine/ 

Water 

120 

240 

7.2 

Mercury 

130 

190 

10.0 

130 

300 

1.5 

190 

560 

0.8 

2 . 2  Droplet  Holography 


Holographic  observations  of  the  droplet  interaction  within  the 
nozzle  contraction  were  produced.  The  holocamera  imaging  optics  were 
configured  with  object-to-hologram  image  magnifications  of  one  and  five 
(Figure  2).  The  large  f ield-of-vlew  optics  (^  »  100  rom  diameter,  lx 
magnification)  were  used  in  the  initial  experiments  in  which  the  droplet 
velocities  were  measured  along  the  nozzle.  Double  exposure  (At=»10gs) 
holography  is  used  to  observe  droplet  velocity.  Smaller  f ield-of-view 
optics  (4>  =  25  mm  diameter,  5x  magnification)  were  used  to  record  high 
resolution  observations  of  the  droplet  breakup  process.  For  this  case, 
external  triggering  was  required  to  center  the  droplet  position  in  the 

hologram.  With  premagnification  (5x)  the  holocamera  spatially  resolves 

•8 

Aura  in  reconstruction.  Since  the  holograms  are  recorded  in  U 
saconds,  droplets  of  order  lOgm  diameter  can  be  resolved  at  velocities 
approaching  10  m/s.  For  these  experiments  the  droplet  velocities  were 
limited  to  about  100  m/s. 

2.3  Laser  Velocimeter 

Droplet  and  gas  phase  velocities  were  obtained  using  laser 
Doppler  veloclraetry  (LDV),  Figure  3.  This  technique  is  based  on  laser 
light  scatter  by  droplets  or  seed  particles  in  the  flow.  A  single  laser 
beam  (Spectra  Physics,  Model  124,  15  ra  watt)  is  split  into  two  equal 
intensity  beams  which  intersect  at  the  desired  measurement  point  in  the 
flow.  Electromagnetic  interference  between  the  two  laser  beams  occurs 
at  the  intersection  region  resulting  in  fringe  formation.  The  fringes 
are  parallel  planes  of  high  Intensity  (constructive  Interference) 
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Figure  3.  Laser  Doppler  Velocimet 


alternating  with  lower  intensity  (destructive  Interference).  For  these 
experiments  the  spacing  between  successive  fringes  was  7pra.  Light 
scattered  by  particles  in  the  flow  produces  signals  modulated  by  the 
fringe  pattern  intensity.  The  electronic  signal  processor  (Macrodyne, 
2000  Series,  LDV  Processor)  measures  the  time  between  successive  fringes 
from  which  particle  velocity  is  obtained.  Distributions  of  a  large 
number  of  events  (typically  10**)  were  obtained  at  each  measurement  point 
in  the  flow  from  which  mean  velocities  were  obtained.  The  droplet  and 
fragment  velocity  distributions  were  measured  with  the  laser  velo- 
cimeter.  The  mean  velocity  is  calculated  at  each  station.  The  droplet 
velocity  profile  was  integrated  to  determine  the  elapsed  time  and  dif~ 
ferentiated  to  determine  the  acceleration.  The  gas/droplet  slip 
velocity  was  used  to  determine  the  Weber  number  (based  on  initial  drop¬ 
let  diameter)  history  along  the  nozzle.  Hence,  both  the  Weber  number  at 
droplet  breakup  and  the  breakup  time  were  determined. 
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3.0 


CONVENTIONAL  LIQUID  EXPERIMENTS 


A  series  of  experiments  was  conducted  to  observe  the  droplet 
deformation  and  fragmentation  processes  with  holography  and  to  determine 
the  droplet  and  fragment  dynamics  with  laser  velocimetry .  The 
experiments  were  conducted  with  three  liquid  types.  For  each  liquid 
with  a  corresponding  surface  tension,  a  parametric  set  of  trajectory 
calculations  were  performed  for  various  nozzle  exit  velocities,  and 
droplet  sizes.  In  an  early  series  of  experiments,  the  rigid  droplet 
trajectory  estimates  were  found  to  grossly  over  estimate  the  Weber 
number  history  when  compared  to  experiment  data.  In  the  experiments, 
the  liquid  droplets  were  found  to  achieve  much  higher  accelerations  Just 
prior  to  breakup.  The  mechanism  for  enhanced  acceleration  was  presumed 
to  be  droplet  flattening,  and  as  a  result  the  droplet  Weber  numbers  were 
lower  than  predicted.  The  gas  velocity  was  increased  so  that  breakup 
occurred.  The  test  conditions  are  listed  in  Table  II;  Including  the 
liquid  type,  droplet  size,  maximum  gas  velocity,  the  slip  velocity  and 
Weber  number  at  breakup,  the  breakup  time,  and  the  period  of  the  first 
natural  frequency  (i.e.  Eqn.  2).  The  breakup  times  were  measured 
relative  to  the  time  at  which  the  Weber  number  exceeded  five. 

3.1  Results 


The  holographic  investigation  revealed  that  the  breakup  mode  for 
nozzle  accelerations  was  similar  to  the  breakup  modes  observed  for  shock 
wave  accelerations.  Initially  the  droplet  flattens  under  the  high 
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pressure  exerted  on  the  stagnation  point  (Figure  4).  As  the  droplet 
flattens,  its  radius  of  curvature  increases  and  the  stagnation  pressure 
is  felt  over  a  larger  area.  The  center  of  the  droplet  is  eventually 
pushed  downstream  of  the  outer  edge  forming  a  thin  membrane  in  the  shape 
of  a  bag  (Figure  5).  When  the  bag  bursts,  very  small  fragments  are 
formed;  however,  considerable  mass  still  remains  in  the  annular  ring. 

The  annular  ring  fails  later  producing  the  larger  fragments  as  has  been 
reported  by  others^^. 

A  glycerine/water  mixture  was  used  to  investigate  viscous 
effects.  The  mixture  was  selected  to  obtain  a  viscosity  10  times  that 
of  water  (Yniixture  ”  cp).  The  significant  effect  of  viscosity  is  to 
modify  the  rate  of  deformation  in  the  streamwlse  and  radial 
directions.  Specifically,  the  rate  of  expansion  in  the  streamwlse  and 
radial  directions  are  observed  to  be  about  equal;  and  as  a  result,  the 
droplets  expand  radially  to  about  7  times  the  original  diameter 
(Figure  6).  For  the  lower  viscosity  case,  radial  expansion  is  limited 
to  about  a  factor  of  four  or  five.  The  mechanism  of  droplet  breakup  is 
not  significantly  altered;  however,  some  details  are  strikingly 
affected.  For  the  higher  viscosity  case,  geometric  symmetry  is 
maintained  throughout  substantial  portions  of  breakup  event;  whereas  for 
the  lower  viscosity  case,  the  deformation  quickly  becomes  asymmetric  or 
random. 

The  laser  veloclmeter  was  used  to  obtain  velocity  distributions 
(Figure  7),  Initially  the  velocity  distribution  is  narrow;  however,  the 
distribution  broader.o  substantially  within  the  breakup  region.  Finally, 
as  a  result  of  the  large  Increase  in  the  mean  velocity  of  the  fragments 
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Figure  4.  Surface  Tension  Effects  on  Droplet  Breakup 
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Figure  5,  Holography  Data,  Alcohol  Droplets  (Series  5,  Table  II) 

(a)  Double  Pulse  (AT  =  25  ys)  Bag  Breakup  Mode 

(b)  Ring  Structure  Observed  Later  in  Breakup  Event 


the  relative  width  of  the  distribution  Is  reduced.  The  broad 
distribution  within  the  breakup  region  should  not  be  confused  with 
turbulence,  but  should  be  associated  with  the  breakup  phenomena.  Fairly 
mild  accelerations  were  observed  until  the  droplet  begins  to  deform  and 
flatten.  The  drag  to  mass  ratio  of  the  flattened  droplet  is 
considerably  higher  than  for  the  spherical  geometry;  and  as  a  result  the 
droplet  experiences  extreme  accelerations  just  before  breakup.  This 
final  acceleration,  of  course,  reduces  the  peak  slip  velocity  and  Weber 
number  and  nust  be  accurately  modeled  if  the  Weber  number  at  breakup  is 
to  be  predicted. 

The  Weber  number  profiles  for  water  droplets  are  shown  in 
Figure  8.  The  position  at  which  breakup  occurs  moves  downstream  if  the 
droplet  size  or  the  nozzle  velocity  is  reduced.  As  a  result  of  reduced 
size  or  velocity  the  peak  Weber  number  is  reduced;  and  eventually,  a 
condition  is  reached  for  which  droplet  breakup  does  not  occur  within  the 
nozzle  contraction.  In  our  experiments  with  water  droplets,  this 
condition  was  reached  only  for  the  smallest  droplet  size,  150  pm,  and 
the  smallest  velocity  rise,  100  m/s,  in  which  case,  the  Weber  number  did 
not  exceed  20  and  failure  did  not  occur  in  the  nozzle  contraction. 

Tlius,  for  a  given  velocity  rise  in  the  nozzle  contraction,  there  exists 
a  boundary  in  the  Weber  number/position  plane  at  which  droplet  breakup 
occurs . 

Since  the  breakup  position  is  dependent  on  the  nozzle  velocity 
rise,  consider  instead  the  time  history  of  loading.  This  is  a  more 
unified  approach  since  it  is  the  loading  duration  which  produces  droplet 
deformations  and  eventually  breakup.  A  zero  time  must  be  selected  and 
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we  have  chosen  che  time  at  which  the  Weber  number  reaches  five.  The 
Weber  number  (Loading)  is  plotted  against  relative  time  and  different 
velocity  rise  and  droplet  size  conditions  compared.  The  relative  time 
has  been  normalized  by  the  period  of  the  first  natural  frequency 
(Eqn.  2). 

The  water  droplet  data  (Figure  9)  fall  along  a  boundary  starting 
at  a  Weber  number  of  about  50  and  time  of  one  and  extending  to  very  low 
Weber  numbers  at  long  times.  One  of  the  alcohol  data  points  is  near  the 
water  boundary,  while  the  others  are  at  much  higher  Weber  numbers  (i.e., 
over  100).  The  glycerine  droplet  data  point  indicates  that  breakup 
occurred  at  a  lower  Weber  number  than  for  water  droplets.  Holographic 
observations  revealed  that  the  glycerine  droplets  flattened  to  a  much 
larger  cross  section  as  a  result  of  the  increased  viscosity.  The  addi¬ 
tional  cross  section  increased  the  drag  and  acceleration  resulting  in  a 
lower  Weber  number  trajectory.  The  fundamental  result  is  Chat  gradual 
loadings  produced  by  nozzle  accelerations  allows  the  droplets  to  reach 
higher  Weber  numbers  than  are  critical  for  impulsive  loading  as  produced 
in  shock  cube  simulations. 


C* 


3.2  Breakup  Parameters 


Weber  number  effect  on  breakup  was  studied  by  changing  both 
droplet  size  and  nozzle  velocity.  The  bag  was  still  observed  Co  expand 
downstream  of  the  ring;  however,  after  the  bag  bursts,  che  ring  exhibits 
more  random  distortion  with  Increasing  Weber  number  (Figure  10), 

Increased  Weber  number  was  observed  to  decrease  the  fragment  size.  The 
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Critical  Weber  Numbers  and  Associated  Time  Durations  of  the  Breakup  Event 
Series  Number  (t.e.  1-8)  refers  to  Table  II 


82-2193 


1. 


LIQUID: 

DIAMETER  (um): 
MAX.  WE.  *: 


t 

tc  =  Pri 

.  ad  ] 


1  =  0?  We  ^=10 


ETHANOL 

ETHANOL 

ETHANOL 

ETHANOL 

164 

265 

164 

265 

47 

38 

90 

98 

1000  >im  BOO  lie*  1000  vm  SOO^jin 


Figure  10.  Weber  Number  Effects  on  Ethanol  Droplet  Breakup 


mean  fragment  size  is  plotted  against  peak  Weber  number  (Figure  11). 

The  mean  fragment  size  based  on  acceleration  (Equation  A)  is  also 
plotted.  The  mean  fragment  size  was  estimated  by  using  a  force  balance 
on  the  fragments.  The  drag  load  was  written  In  terms  of  the  slip 
velocity,  fragment  size,  and  a  drag  coefficient  based  on  the  Reynolds 
number.  The  mass  was  written  in  terms  of  the  fragment  size  and  density 
(the  fragments  are  assumed  spherical,  a  fact  confirmed  by  holographic 
data);  and  the  acceleration  was  measured. 

The  force  balance  resulted  in  an  equation  for  the  fragment  size, 
1^,  as  follows: 

F  =  MA 
or 

n  ^  IT  ^ 

■Pgi'  A 

where 

q  =  dynamic  pressure  based  on  the  slip  velocity 
Cq  =  drag  coefficient 
p  -  liquid  density 
A  =•  acceleration 


leaving 


(A) 
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Symbols  =  LDV  Data;  Closed  Symbols  =  Holography  Data 


For  the  alcohol  droplets  both  estimates  of  the  droplet  size  were  In 
agreement.  The  fragment  size  for  water  and  glycerine/water  mixtures  are 
also  plotted;  however  no  imaging  data  was  available  for  comparison. 

The  large  fragment  size  estimated  for  water  droplets  was  not  expected 
and  may  indicate  that  perhaps  Che  larger  fragments  were  not  spherical. 
The  small  fragment  size  for  glycerine/water  mixture  may  be  a  result  of 
the  initial  droplet  expansion  to  about  seven  times  original  diameter 
and,  therefore,  the  bag  membrane  thickness  must  be  reduced.  Hence,  the 
effect  of  Increased  liquid  viscosity  was  reduced  fragment  size. 
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4.0  LIQUID  METAL  EXPERIMENTS 


A  series  of  experiments  were  designed  to  investigate  the  breakup 
dynamics  of  high  surface  tension  liquid  metals.  Many  liquid  metals  were 
identified,  but  mercury  was  selected  as  it  was  compatible  with  the 
droplet  generation  technique,  and  with  room  temperature  air  experi¬ 
ments.  The  aerodynamics  of  the  experiment  were  based  on  trajectory 
estimates  for  mercury  droplets  from  100  to  600  ym.  Since  the  liquid-to¬ 
gas  density  ratio  was  so  large,  the  droplet  velocities  were  assumed  to 
be  small  and  the  slip  velocity  was  taken  as  the  gas  velocity.  This  fact 
was  verified  in  the  experiments  when  the  droplet  velocity  was  limited  to 
less  than  10  percent  of  the  gas  velocity.  For  breakup  (We  ■  30.),  the 
slip  velocity  or  gas  velocity  requirement  ranged  from  270  m/s  to  109  m/s 
for  droplet  diameters  of  100  ym  to  600  ym,  respectively.  A  series  of 
experiments  were  defined  in  which  the  nozzle  would  be  operated  at  sonic 
velocity  and  the  droplet  size  would  be  varied.  The  2-dlmenslonal  nozzle 
exhaust  duct  was  modified  to  achieve  sonic  velocity  in  the  throat.  The 

nozzle  would  start  and  was  operated  at  a  stagnation  pressure  of  1.84  x 
5  2 

10  N/M  .  A  droplet  generator  with  a  130  ym  diameter  orifice  was 
operated  at  10  kHz  and  1.5  kHz,  to  produce  190  ym  and  300  ym  diameter 
droplets,  respectively;  and  a  second  orifice  size  ((^  ■  190  ym)  was 
operated  at  840  Hz  to  produce  560  ym  diameter  droplets.  The  test  matrix 
is  shown  in  Table  II,  The  pulsed  laser  holography  was  used  to  record 
droplet  Images  in  various  states  along  the  nozzle. 
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4.1 


Results 


The  droplet  breakup  process  was  very  similar  for  the  three 
droplet  sizes,  except  that  the  breakup  moved  to  a  higher  gas  velocity 
position  in  the  nozzle  for  smaller  droplets.  The  breakup  mechanism 
(Figure  12)  is  observed  to  be  the  Umbrella  Mode  characteristic  of  the 
Weber  number  range  of  70-100  for  conventional  liquids.  In  this  breakup 
mode,  the  droplet  initially  flattens  and  then  expands  to  larger 
dimension  (3-5X)  forming  the  stem  in  the  center  and  a  sheet  between  the 
stem  and  the  outer  ring.  The  sheet  bursts  quickly  constituting  the 
primary  breakup  event.  The  outer  ring  fragments  soon  after  the  sheet 
bursts.  The  stem  is  stable  at  this  point  and  remains  secure  until  a 
higher  slip  velocity  is  reached.  The  sheet  and  ring  fragments  are  much 
smaller  and  accelerate  away  from  the  massive  stem  (Figure  13).  The  stem 
is  rather  cylindrical  and  aligned  with  the  flow.  The  frontal  diameter 
of  the  stem  is  about  30-40  percent  of  the  original  droplet  size  and 
deforms  as  higher  velocities  are  reached  in  the  nozzle.  Secondary 
breakup  (Figure  12)  is  observed  as  the  stem  approaches  the  throat.  The 
primary  fragments  travel  downstream  of  the  stem  during  the  secondary 
breakup  cycle,  and  at  the  throat  the  most  distant  fragments  (which  are 
also  the  smallest)  are  about  25.  mm  downstream  of  the  stem.  Finally, 
after  the  stem  has  failed,  the  breakup  is  essentially  complete  and  most 
of  the  smaller  fragments  are  spherical.  The  larger  fragments  arise  from 
the  stem  and  ring  and  the  smaller  fragments  arise  from  the  sheet. 

The  velocity  and  Weber  number  profiles  (Figure  14)  depict  the 
velocities  of  the  initial  droplet,  the  primary  fragments,  and  the 
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Figure  15.  Primary  Fragments  Resulting  from  Metal  Droplet  Breakup 
(Series  11,  Table  II) 


Series 


stem.  The  position  (upstream  of  the  throat)  for  the  primary  breakup 
moves  from  8.0  mm  to  14.0  mm  for  the  initial  droplet  sizes,  190  pm  to 
560  uni,  respectively.  The  Weber  number  is  based  on  the  slip  velocity 
and  initial  droplet  size;  then,  for  the  stem,  the  Weber  number  is  based 
on  an  approximate  frontal  diameter  (typically  =*  .30  The  initial 

droplets  all  fail  in  the  Weber  number  range  of  15-20,  and  the  stems  fall 
in  a  similar  range.  The  Weber  number  histories  (Figure  15)  of  all  three 
sizes  overlap  to  within  the  accuracy  of  the  data,  when  the  time  is 
normalized  by  the  period  of  the  first  natural  frequency  (Eqn.  2).  This 
period  for  acoustic  oscillation  and  the  breakup  time  are  essentially 
equivalent . 

The  fragment  size  distribution  was  determined  from  high  magnifi¬ 
cation  photographs  of  the  reconstructed  holographic  Images  (ten  events 
were  sampled.).  An  estimate  of  the  equivalent  spherical  diameter  was 
used  for  nonspherlcal  fragments.  The  fragment  number  density 
distribution  (Figure  16)  ranges  from  5-50  u®  and  the  sampling  volume  for 
each  droplet  size  bin  is  corrected  for  the  depth  of  field.  The  mean 
fragment  size  is  about  10  percent  of  the  original  droplet  diameter. 

Each  fragment  number  density/size  distribution  is  a  result  of  a  single 
droplet  breakup  event  in  the  nozzle  contraction.  The  sample  is  obtained 
in  an  instant  In  time,  but  is  spatially  averaged  along  the  nozzle. 
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SUMMARY 


Droplet  breakup  experiments  have  been  performed  in  laboratory 
scale  nozzle  contraction  (throat  size  =  12  mm  square).  Room  temperature 
air  was  used  for  the  test  gas  and  a  conventional  liquid  and  a  liquid 
metal  was  used  to  produce  the  droplets.  The  droplet  size  was  varied 
from  100  um  to  600  pm  and  the  nozzle  velocity  rise  was  varied  from 
100  m/s  to  300  m/s.  Droplet  velocity  profiles  along  the  nozzle  were 
resolved  with  laser  velocimetry  and  with  double  pulsed  laser  holography 
which  also  provided  graphic  images  of  the  droplet  breakup  mechanisms  and 
of  the  droplet  fragments.  The  droplet  velocity  profiles  were  used  to 
determine  the  Weber  number  profiles  and  histories.  Thus,  the  critical 
Weber  number  and  a  measure  of  the  duration  of  the  droplet  breakup  event 
was  determined. 

The  breakup  mechanism  for  conventional  liquid  droplets  was,  in 
general,  equivalent  to  the  Bag  type  previously  observed  in  shock  tube 
experiments.  The  critical  Weber  number  and  breakup  time  data  revealed 
that  when  breakup  occurs  in  a  short  time  that  higher  V/eber  numbers  are 
required;  and  conversely  when  breakup  occurs  in  a  long  time  that  lower 
Weber  numbers  are  required.  This  trend  was  observed  for  both  water  and 
alcohol  droplets.  For  water  droplets  we  observed  a  critical  Weber 
number  of  about  50  for  breakup  times  of  about  2-4  which  then  decayed  to 
a  data  point  for  which  breakup  did  not  occur  in  the  nozzle  contrac¬ 
tion.  One  of  the  alcohol  data  points  fell  on  the  water  boundary,  while 
the  others  reached  the  highest  critical  Weber  numbers  observed  in  these 
experiments.  The  data  obtained  for  the  more  viscous  liquid 


-35- 


84-2193-11/D9 


(a  glycerine/water  mixture)  revealed  the  effect  of  increased  viscosity 
to  be  increased  acceleration.  Viscous  liquid  droplets  were  observed  to 
flatten  to  almost  seven  times  their  original  diameter  before  breakup 
occurred.  The  Increased  acceleration  resulting  from  the  larger  cross 
section  reduces  the  slip  velocity  or  Weber  number  achieved  within  the 
nozzle.  An  additional  effect  of  Increased  viscosity  is  a  reduction  in 
the  fragment  size.  The  ratio  of  the  mean  fragment  to  initial  droplet 
size  was  determined  for  various  Weber  number  failures.  For  alcohol  and 
water  droplets,  the  ratio  varied  from  40%  at  low  Weber  number,  50,  to 
about  10%  at  high  Weber  number,  120.  Interesting,  the  high  viscosity 
data  point  reveals  a  fairly  small  fragment  size  for  the  low  Weber 
number,  40. 

The  liquid  metal  droplet  experiments  required  higher  dynamic 
pressures  for  breakup  and  as  a  result  the  nozzle  was  operated  at  the 
sonic  velocity.  The  breakup  mechanism  was  observed  to  be  the  "Umbrella" 
type.  This  mechanism  has  been  observed  for  mercury  droplets  In  shock 
tube  experiments,  but  only  for  Weber  numbers  over  80.  The  breakup  mech¬ 
anism  for  lower  Weber  number  shock  cube  experiments  "Bag"  type.  A  large 
stem  resulted  from  all  primary  breakup  events  which  eventually  failed 
further  down  the  nozzle  as  higher  slip  velocities  are  achieved. 

Secondary  breakup  is  an  important  feature  of  the  metal  droplet  breakup 
mechanism.  In  contrast  to  the  conventional  liquid  results,  the  Weber 
number  histories  for  liquid  metals  collapsed  to  a  single  curve 
terminated  by  the  primary  breakup  event.  All  three  droplet  sizes 
experienced  breakup  at  the  same  Weber  number,  20.  Also,  the  duration  of 
the  breakup  event  was  observed  to  be  a  constant  when  normalized  by  the 
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natural  period.  The  fragment  size  distribution  was  observed  to  have  a 
mean  value  of  about  ten  percent  the  original  droplet  size.  Since  the 
critical  Weber  number  did  not  vary  one  might  expect  the  relative  mean 
fragment  size  not  to  vary. 
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AbSi-ract 


The  breakup  of  AL/AL20-J  agglomerates  in  solid 
propellant  rocket  nozzles  is  a  critical  process 
effecting  combustion  efficiency  and  two-phase  flow 
losses.  Droplet  breakup  experiments  in  aerody¬ 
namic  nozzle  contractions  have  been  conducted  in 
which  conventional  liquids  and  higher  surface  ten¬ 
sion  liquids  (mercury)  were  examined.  Pulsed 
laser  holography  has  provided  droplet  and  fragment 
observations  with  resolution  not  previously 
obtained.  Laser  veloclmetry  has  provided  droplet 
dynamics  data  revealing  dramatic  accelerations 
prior  to  breakup.  Conventional  liquid  droplets 
were  observed  to  survive  to  much  higher  Weber  num¬ 
bers  (of  order  100)  than  those  considered  to  be 
critli'ul  for  breakup.  In  contrast,  the  liquid 
metal  droplets  were  observed  to  fail  at  fixed 
Weber  number  of  about  15-20  and  Che  duration  of 
the  failure  event  was  found  to  scale  with  a  time 
scale  based  on  the  liquid  properties. 

Introduction 


Combustion  efficiency  of  aluminized  propel¬ 
lants  in  solid  rocket  motors  is  reduced  by  incom¬ 
plete  aluminum  combustion  and  two-phase  nozzle 
flow  losses.  The  combustion  of  aluminized  propel¬ 
lants  can  produce  large  AI/AI2OJ  agglomerates. 
The  agglomerate  dynamics  within  the  combustion 
chamber  and  nozzle  have  a  significant  influence  on 
the  overall  combustion  efficiency  of  the  motor**'. 
Agglomerates  are  subjected  to  large  aerodynamic 
loads  within  the  rocket  nozzle  where  the  gas  phase 
experiences  a  more  rapid  acceleration  than  the 
agglomerates.  The  drag  load  deforms  the  agglomer¬ 
ates  and  if  cf  sufficient  magnitude  result  in 
breakup.  The  smaller  fragments  have  faster  velo¬ 
city  and  thermal  equilibrium  times  and  have  higher 
combustion  rates.  For  maximum  combustion  effi¬ 
ciency  the  aluminum  fragments  must  completely  oxi¬ 
dize,  and  achieve  kinetic  and  thermal  equilibrium 
with  the  gas  phase.  As  a  direct  result  of  agglom¬ 
erate  breakup,  Che  aluminum  combustion  rate  is 
Increased,  and  the  thermal  energy  released  is  more 
efficiently  transferred  into  exhaust  kinetic 
energy. 

Photographic  observations  obtained  in  windowed 
rocket  motors”  and  combustion  bombs””***  indicate 
relatively  large  agglomerates  (100pm  -  500pm)  are 
formed  on  the  propellant  surface  and  entrained  in 
Che  combustion  flow;  however,  particle  size 
measureTCDts  obtained  from  sampling  the  exhaust 
plumes*  indicate  small  mean  particle  diameters 

(<10pm).  These  small  exhaust  plume  particles 
apparently  result  from  the  breakup  of  the  larger 
agglomerates  during  the  nozzle  expansion  process. 
Observations  of  agglomerate  breakup  in  a  labora¬ 
tory  scale  rocket  nozzle  revealed  an  adequate  cor¬ 
relation  with  Weber  number;  however,  neither  the 


physical  process  of  broakiyi  nor  the  fragment  size 
distribution  was  resolved*  . 

This  research  has  the  objective  of  obtaining 
physical  data  to  characterize  the  mechanisms  of 
aerodynamic  droplet  breakup.  An  experiment  has 
been  completed  in  which  conventional  liquids 
(water,  alcohol  and  glycerlne/water  itlxtures)  and 
a  liquid  metal  (mercury)  was  studied.  The  primary 
goal  of  Che  conventional  liquid  experiments  was  to 
examine  the  effect  of  liquid  properties  (viscosity 
and  surface  tension)  on  the  breakup  rn'chanlsm, 
time  scale,  and  fragment  size  distribution.  The 
goal  of  the  liquid  metal  experlmercs  was  to  exam¬ 
ine  the  effect  of  the  much  higher  surface  tension 
more  characteristic  of  liquid  aluminum. 

A  key  element  of  the  experimental  effort  is 
Che  use  of  nonintruslve  laser  diagnostics  includ¬ 
ing  pulsed  laser  holography  and  laser  Doppler 
veloclmetry.  The  exceptional  tempcral  and  spatial 
resolution  of  pulsed  laser  holography  provided  the 
ability  to  resolve  the  mechanism  of  breakup  and 
the  size  distribution  of  the  fragments.  Laser 
Doppler  veloclmetry  was  used  to  determ. ne  drop 
velocity  distributions  along  the  nozzle  revealing 
the  rapid  acceleration  of  the  flattenej  droplets 
and  then  surprisingly  the  milder  acceleration  of 
the  fragments.  The  high  drag  to  macs  ratio  of  the 
flattened  droplets  was  expected  and  re8o)ved  dur¬ 
ing  the  initial  phase  of  the  breakup  process. 

Experiment  Considerations 

Aerodynamic  droplet  breakup  is  characterized 
by  Che  magnitude  and  duration  of  the  aerodynamic 
forces.  For  nozzle  contractions,  the  droplets  are 
loaded  as  a  result  of  the  higher  gas  phase  accel¬ 
eration.  The  droplets  respond  to  the  aerodynamic 
forces  by  deforroing  and  accelerating.  The  initial 
acceleration  scales  with  the  initial  diameter; 
however,  as  the  pressure  forces  flatten  the  drop¬ 
let,  increasing  its  cross-sectional  area,  the  drag 
load  and  the  acceleiatlon  Increases  rapidly.  This 
sequence  of  events  can,  of  course,  lead  to  catas¬ 
trophic  deformation  and  droplet  breakup  if  the 
slip  velocity  is  of  sufficient  magnitude  and  dura¬ 
tion. 


The  liquid  surface  tension  Is  used  to  scale 
the  aerodynamic  load  forming  the  Weber  number: 


We 
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The  duration  of  the  aerodynamic  loading  can  be 
scaled  by  the  natural  period  for  hydrodynainlc 
oscillation,  which  Is  a  measure  of  the  droplet 
response  time: 


For  nozzle  accelerations,  an  additional  time  scale 
exists;  namely  the  tine  required  for  velocity 
equilibrium,  T  .  Assuming  the  droplet  Reynolds 
number  relative  to  the  gas  phase  la  sufficiently 
high  (i.e..  Re  >  10^  ),  the  drag  coefficient  for 
rigid  droplets  Is  approximately  unity  and  the  cor¬ 
responding  time  scale  la 


However,  since  Che  liquid  droplets  are  not  rigid, 
the  droplet  cross  section  and  the  associated  drag 
load  Increase  rapidly  within  the  nozzle  contrac¬ 
tion.  Therefore,  the  time  scale  for  liquid  drop¬ 
let  velocity  equilibrium  Is  no  doubt  faster  than 
Che  rigid  droplet  estlaiate. 

A  few  liquid  types  were  selected  to  examine 
Che  effect  of  liquid  properties  such  as  surface 
tension  and  viscosity.  For  the  conventional 

liquids,  water  and  alcohol  were  selected  for  Che 
surface  tension  series;  and  water  and  a  glycerine 
water  (57/43)  mixture  were  selected  for  the  vis¬ 
cosity  aeries.  The  mixture  was  selected  Co 
achieve  a  viscosity  ten  tines  that  of  water, 

A  second  series  of  experiments  were  designed 
CO  exsoUne  the  breakup  dynamics  of  high  surface 
tension  liquid  metals.  Mercury  was  selected 

because  of  its  compatibility  with  the  droplet  gen¬ 
eration  technique  and  with  room  temperature  sir 
flow. 

Nozzle  Facility 

An  experiment  was  designed  In  which  droplet 
breakup  could  be  Investigated  in  a  nozzle  contrsc- 
Cion.  The  nozzle  wss  designed  to  allow  droplet 
injection  upstream  of  the  contraction,  Windows 
formed  the  sides  of  the  two-dimensional  nozzle 
allowing  observation  of  the  droplet  Interaction, 
The  gas  was  delivered  from  a  settling  chamber 
through  a  sintered  metal  throttle  to  the  nozzle 
plenum  where  the  droplets  were  Injected,  The  noz¬ 
zle  area  ratio  was  eight.  Fig,  1,  For  the  conven¬ 
tional  liquid  experiments,  the  nozzle  was  operated 
at  subsonic  exit  velocities.  However,  for  the 
liquid  metal  experiments,  high  dynamic  pressures 
were  required  to  achieve  droplet  breakup.  Higher 
droplet  loading  was  achieved  by  Increasing  the  gas 
velocity  and  hence  the  slip  velocity.  In  fact, 
the  nozzle  exhaust  was  modified  to  allow  sonic 
flow  In  the  nozzle  throat  for  the  liquid  metal 
experiments,  although  droplet  breakup  occurred  up¬ 
stream  of  the  throat. 

Piezoelectric  droplet  generators  were  con¬ 
structed  to  produce  highly  monodlsperse  droplet 
streams.  The  basic  design  wss  that  of  modern  Ink¬ 
jet  printers.  The  complete  droplet  generator  Is 
composed  of  a  liquid  supply  tank,  an  electrical 
driver,  and  the  Ink-jet  devices  (Figure  1).  The 


size  of  the  droplets  generated  by  our  oevlces  was 
variable  from  about  1.5  to  3.0  times  the  orifice 
diameter  by  decreasing  the  excitation  frequency. 
For  conventional  liquids,  droplets  were  generated 
from  160  pm  to  400  pm  in  diameter,  and  for  mer¬ 
cury,  droplets  were  generated  from  190  pm  to 
560  pm.  The  droplet  generator  operating  condi¬ 
tions  are  listed  In  Table  1  for  the  conventional 
liquids  and  mercury, 

TABLE  I,  DROPLET  GENERATOR  OPERATING  CONDITIONS 


LIQUID 

DIAMETER 

EXCITATION 

FREQUENCY 

(kHz) 

TIP 

(iim) 

DROPLET 

(um) 

Water 

lOU 

160 

10.6 

120 

220 

8.0 

190 

410 

0.8 

Alcohol 

lOU 

164 

11.2 

150 

265 

9.0 

Glycerine/ 

120 

240 

7.2 

Hater 

Mercury 

130 

190 

IC.O 

130 

300 

l.S 

190 

560 

0.8 

Droplet  Holography 


Holographic  obaervatlona  of  the  droplet 
interaction  within  the  nozzle  contraction  were 
produced.  The  holocaoera  imaging  optics  were 
configured  with  object-co-hologram  image  oagnl- 
flcations  of  one  and  five.  The  large  fleld-of- 
view  optics  (^  ■  100  OD  diameter,  lx  magnifica¬ 
tion)  were  used  in  the  initial  experiments  In 
which  the  droplet  velocities  were  measured  along 
the  nozzle.  Double  exposure  (AtslOpa)  holography 
is  used  to  observe  droplet  velocity,  Soialler 
f Icld-of-vleu  optics  (^  ■  25  mo  diameter,  5x 
magnification)  were  used  to  record  high  resolution 
observations  of  the  droplet  breakup  process.  For 
this  case,  external  triggering  was  required  to 
center  the  droplet  position  In  the  hologram.  With 
premsgnlf Icatlon  (5x)  the  holocamera  spatially 
resolves  4pm  In  reconstruction.  Since  the  holo¬ 
grams  are  recorded  In  10  seconds,  droplets  of 
order  10pm  dlswter  can  be  resolved  at  velocities 
approaching  10''  m/s.  For  these  experiments  the 
droplet  velocities  were  Hod. ted  to  about  100  o/s. 

Laser  Veloclmeter 

Droplet  and  gas  phase  velocities  were  obtained 
using  laser  Doppler  veloclmetry  (LDV),  This  tech¬ 
nique  Is  based  on  laser  light  scatter  by  droplets 
or  seed  particles  in  the  flow.  A  single  laser 
beam  (Spectra  Physics,  Model  124,  15  m  watt)  la 
spilt  Into  two  equal  Intensity  beams  which  Inter¬ 
sect  at  the  desired  measurement  point  In  the  flow. 
Electroiaagnetlc  Interference  between  the  two  laser 
beams  occurs  at  the  Intersection  region  lesulclng 
In  fringe  formation.  The  fringes  are  parallel 
planes  of  high  intensity  (constructive  Interfer¬ 
ence)  alternating  with  lower  Intensity 
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(descrurtlve  Interference).  For  these  experiments 
the  spacing  between  successive  fringes  was  7uni, 
Light  scattered  by  particles  In  the  flow  produces 
signals  modulated  by  the  fringe  pattern  Intensity. 
The  electronic  signal  processor  (Macrodyne,  2000 
Series,  LDV  Processor)  measures  the  time  between 
successive  fringes  from  which  particle  velocity  Is 
obtained.  Distributions  of  a  large  number  of 
events  (typically  10'*  )  were  obtained  at  each 
measurement  point  In  the  flow  from  which  mean 
velocities  were  obtained.  The  droplet  snd  fragment 
velocity  distributions  were  measured  with  the  laser 
veloclmeter.  The  mean  velocity  Is  calculated  at 
each  station.  The  droplet  velocity  profile  was 
Integrated  to  determine  the  elapsed  time  and  dif¬ 
ferentiated  to  determine  Che  acceleration.  The 
gas/droplet  slip  velocity  was  used  to  determine  the 
Weber  number  (based  on  initial  droplet  diameter) 
history  along  the  nozzle.  Hence,  both  the  Weber 
number  at  droplet  breakup  and  the  breakup  time  were 
determined. 


Conventional  Liquid  Experiments 


A  series  of  experiments  were  conducted  to 
observe  the  droplet  deformation  and  fragmentation 
processes  with  holography  and  to  determine  the 
droplet  and  fragment  dynamics  with  laser  velotlm- 
etry.  The  experiments  were  conducted  with  three 
liquid  types.  For  each  liquid  with  a  corresoondlng 
surface  tension,  a  parametric  set  of  trajectory 
calculations  were  performed  for  various  nozzle  exit 
velocities,  and  droplet  sizes.  In  an  early  series 
of  experiments,  the  rigid  droplet  trajectory  esti¬ 
mates  were  found  to  grossly  over  estimate  the  Weber 
number  history  when  compared  to  experiment  data. 
In  the  experiments,  the  liquid  droplets  were  found 
to  achieve  much  higher  accelerations  Just  prior  to 
breakup.  The  mechanism  for  enhanced  acceleration 
was  presumed  to  be  droplet  flattening,  and  as  a 
result  the  droplet  Weber  numbers  were  lower  than 
predicted.  The  gas  velocity  was  Increased  so  that 
breakup  occurred.  The  test  conditions  are  listed 
In  Table  II;  Including  the  liquid  type,  droplet 


size,  maximum  gas  velocity,  the  slip  velocity  and 
Weber  number  at  breakup,  the  breakup  time,  and  the 
period  of  the  first  natural  frequency  (l.e. 
Eqn.  2).  The  breakup  times  were  measured  relative 
to  the  time  at  which  the  Weber  number  exceeded 
five. 

Results  ' 

The  holographic  Investigation  revealed  that  the 
breakup  mode  for  nozzle  accelerations  was  similar 
to  the  breakup  modes  observed  for  shock  wave  accel¬ 
erations.  Initially  the  droplet  flattens  under  Che 
high  pressure  exerted  on  the  stagnation  point.  As 
the  droplet  flattens.  Its  radius  of  curvature 
Increases  and  the  scagnoclon  pressure  Is  felt  over 
a  larger  area.  The  center  of  the  droplet  Is  event¬ 
ually  pushed  downstream  of  the  outer  edge  forming  a 
thin  membrane  In  the  shape  of  a  bag  (Figure  2). 
When  the  bag  bursts,  very  small  fragments  are 
formed;  however,  considerable  mass  still  remains  In 
Che  annular  ring.  The  annular  ring  falls  later 
producing  the  larger  fragments  as  has  been  reported 
by  others*  . 

The  laser  veloclmeter  was  used  to  obtain  velo¬ 
city  distributions  (Figure  3).  Initially  Che 

velocity  distribution  Is  narrow;  however,  rhe  dis¬ 
tribution  broadens  substantially  within  the  breakup 
region.  Finally,  as  a  result  of  the  large  Increase 
in  the  mean  velocity  of  the  fragments  the  relative 
width  of  the  distribution  Is  reduced.  The  broad 
distribution  within  the  breakup  region  should  not 
be  confused  with  turbulence,  but  should  be  associ¬ 
ated  with  the  breakup  phenomena.  Fairly  mild 
accelerations  were  observed  until  the  droplet 
begins  to  deform  and  flatten.  The  drag  to  mass 
ratio  of  the  flattened  droplet  is  considerably 
higher  than  for  the  spherical  geometry;  and  as  a 
result  the  droplet  experiences  extreme  accelera¬ 
tions  Just  before  breakup.  This  final  accelera¬ 
tion,  of  course,  reduces  the  peak  slip  velocity  and 
Weber  number  and  must  be  accurately  modeled  If  the 
Weber  number  at  breakup  Is  to  be  predicted. 


3 


,  w  ^  v.' v'T'v  t  ■ 


The  Weber  number  profiles  for  water  droplets 
are  shown  in  Figure  4,  The  position  at  which 
breakup  occurs  moves  downstream  If  the  droplet  size 
or  the  nozzle  velocity  Is  reduced^  As  a  result  of 
reduced  size  or  velocity  the  peek  Weber  number  is 
reduced;  and  eventually,  a  condition  Is  reached  for 
which  droplet  breakup  does  not  occur  within  the 
nozzle  contraction.  In  our  experiments  with  water 
droplets,  this  condition  was  reached  only  for  the 
smallest  droplet  size,  ISO  go,  and  the  smallest 
velocity  rise,  100  m/s.  In  which  case,  the  Weber 
number  did  not  exceed  20  and  failure  did  not  occur 
In  the  nozzle  contraction.  Thus,  for  a  given 
velocity  rise  In  the  nozzle  contraction,  there 
exists  a  boundary  In  the  Weber  number/posltlon 
plane  at  which  droplet  breakup  occurs. 

Since  the  breakup  position  Is  dependent  on  the 
nozzle  velocity  rise,  consider  Instead  the  time 
history  of  loading.  This  Is  a  more  unified 
approach  since  It  is  the  loading  duration  which 
produces  droplet  deformations  and  eventually  break¬ 
up.  A  zero  time  must  be  selected  and  we  have 
chosen  the  time  at  which  the  Weber  number  reaches 
five.  The  Weber  number  (loading)  Is  plotted 
against  relative  time  and  different  velocity  rise 
and  droplet  size  conditions  compared.  The  relative 
time  has  been  normalized  by  the  period  of  the  first 
natural  frequency  (Eqn.  2). 

The  water  droplet  data  (Figure  S)  fall  along  a 
boundary  starting  at  a  Weber  number  of  about  50  nd 
time  of  one  and  extending  to  very  low  Weber  numbers 
at  long  times.  One  of  the  alcohol  data  points  Is 
near  the  water  boundary,  while  the  others  are  at 
much  higher  Weber  numbers  (l.e,,  over  100).  The 
glycerine  droplet  data  point  Indicates  that  breakup 
occurred  at  a  lower  Weber  number  than  for  water 
droplets.  Holographic  observations  revealed  that 
the  glycerine  droplets  flattened  to  a  imch  larger 
cross  section  as  a  result  of  the  increased  vis¬ 
cosity.  The  additional  cross  section  Increased  the 
drag  and  acceleration  resulting  in  a  lower  Weber 
number  trajectory.  The  fundamental  result  Is  that 
gradual  loadings  produced  by  nozzle  accelerations 
allows  the  droplets  to  reach  higher  Weber  numbers 
than  are  critical  for  Impulsive  loading  as  produced 
In  shock  tube  simulations. 

Breakup  Parameters 

Weber  number  effect  on  breakup  was  studied  by 
changing  both  droplet  size  and  nozzle  velocity. 
The  bag  was  still  observed  to  expand  downstream  of 
the  ring;  however,  after  the  bag  bursts,  the  ring 
exhibits  more  random  distortion  with  Increasing 
Weber  number.  Increased  Weber  number  was  observed 
to  decrease  the  fragment  size.  The  mean  fragment 
size  Is  plotted  against  peak  Weber  number  (Fig, 6). 
The  mean  fragment  size  based  on  acceleration  (Equa¬ 
tion  4)  is  also  plotted.  The  mean  fragment  size 
was  estimated  by  using  a  force  balance  on  the  frag¬ 
ments.  The  drag  load  was  written  in  terms  of  the 
slip  velocity,  fragment  size,  and  a  drag  coeffi¬ 
cient  based  on  the  Reynolds  number.  The  mass  was 
written  In  terms  of  the  fragment  size  and  density 
(the  fragments  are  assumed  spherical,  a  fact  con¬ 
firmed  by  holographic  data);  and  the  acceleration 
was  measured. 
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q  ■  dynamic  pressure  based  on  the  slip  velocity 
Cjj  ”  drag  coetficient 
p  -  liquid  density 
A  “  acceleration 
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For  the  alcohol  droplets  both  estimates  of  the 
droplet  size  were  in  agreement.  The  fragment  size 
for  water  and  glycerine/water  mixtures  are  also 
plotted;  however  no  imaging  data  was  available  for 
comparison.  The  large  fragment  size  estimated  for 
water  droplets  wss  nor  expected  and  may  Indicate 
that  perhaps  the  larger  fragments  were  not  spher¬ 
ical.  The  small  fragment  size  for  glycerine/water 
mixture  may  be  a  result  of  the  Initial  droplet 
expansion  to  sbout  seven  times  original  diameter 
and,  therefore,  the  bag  membrane  thickness  ousc  be 
reduced.  Hence,  the  effect  of  Inci sased  liquid 
viscosity  was  reduced  fragment  size. 

Liquid  Metal  Experiments 

A  series  of  experiiaents  were  designed  to 
Investigate  the  breakup  dynamics  of  high  surface 
tension  liquid  metals.  Many  liquid  metals  were 
Identified,  but  mercury  was  selected  as  It  was  com¬ 
patible  with  the  droplet  generation  technique,  and 
with  room  temperature  air  experlmenta.  The  aerody¬ 
namics  of  Che  experiment  were  based  on  trajectory 
estimates  for  mercury  droplets  from  100  to  600  gm. 
Since  the  llquld-co-gss  density  ratio  was  so  large, 
Che  droplet  velocities  were  assumed  to  be  small  and 
Che  slip  velocity  was  taken  as  the  gas  velocity. 
This  fact  was  verified  In  the  experlmenta  vrhen  the 
droplet  velocity  waa  limited  to  less  than  10  per¬ 
cent  of  the  gas  velocity.  For  breakup  (We  “  30.), 
the  slip  velocity  or  gas  velocity  requirement 
ranged  from  270  m/a  to  109  m/s  for  droplet  diame¬ 
ters  of  100  gm  CO  600  gm,  respectively.  A  series 
of  experiments  were  defined  In  which  Che  nozzle 
would  hr  operated  at  sonic  velocity  and  the  droplet 
size  would  be  varied.  The  2-dlmen8lonal  nozzle 
exhaust  duct  was  modified  to  achieve  sonic  velocity 
In  the  throat.  The  nozzle  would  start  and  was 
operated  at  a  stagnation  pressure  of  1.84  x  10^ 
N/M*.  A  droplet  generator  with  a  130  gm  diameter 
orifice  was  operated  at  10  kHz  and  1.5  kHz,  to  pro¬ 
duce  190  gm  and  300  gm  diameter  droplets,  respec¬ 
tively;  and  a  second  orifice  size  (^  •  190  gm)  was 
operated  at  840  Hz  to  produce  560  gm  diameter  drop¬ 
lets.  The  test  matrix  la  shown  In  Table  11,  The 
pulsed  laser  holography  was  used  to  record  droplet 
images  In  various  states  along  the  nozzle. 


The  force  balance  resulted  in  an  equation  for 
the  fragment  size,  as  follows: 


The  droplet  breakup  process  was  very  similar 
for  the  r.'iree  droplet  sizes,  except  that  the  break¬ 
up  moved  to  a  higher  gas  velocity  position  in  the 
nozzle  for  smaller  droplets.  The  breakup  mechanism 


(Figure  7)  la  observed  Co  be  Che  Umbrella  Mode 
characcerisclc  of  Che  Weber  number  range  of  70—100 
for  convencional  liquids*  In  Chls  breakup  mode* 
Che  droplet  iniiially  flaccens  and  Chen  expands  Co 
larger  dimension  (3-3X)  forming  Che  sCem  in  Che 
center  and  a  sheet  between  the  seem  and  Che  outer 
ring.  The  sheet  bursts  quickly  constituting  Che 
pripsry  breakup  event.  The  outer  ring  fragments 
soon  after  the  sheet  bursts.  The  stem  is  stable  at 
this  point  and  remains  secure  until  a  higher  slip 
velocity  is  reached.  The  sheet  and  ring  fragments 
are  much  suller  and  accelerate  away  from  Che 
massive  stem  (Figure  8).  The  stem  is  rather  cylin¬ 
drical  and  aligned  with  the  flow.  The  frontal 
diameter  of  Che  stem  is  about  30-40  percent  of  the 
original  droplet  size  and  defr.ms  as  higher  veloci¬ 
ties  are  reached  in  Che  nozzle.  Secondary  breakup 
(Figure  7)  is  observed  as  the  stem  approaches  Che 
throat.  The  primary  fragments  travel  downstream  of 
Che  stem  during  Che  secondary  breakup  cycle,  and  at 
Che  throat  the  most  distant  fragments  (which  are 
also  the  smallest)  are  about  2S.  am  downstream  of 
the  stem.  'laally,  after  the  stem  has  failed,  the 
breakup  is  essentially  complete  and  most  of  the 
smaller  fragments  are  spherical.  The  larger  frag¬ 
ments  arise  from  Che  stem  and  ring  and  Che  smaller 
fragments  arise  from  the  sheet. 

The  velocity  and  Weber  number  profiles  (Fig¬ 
ure  9)  depict  the  velocities  of  Che  initial  drop¬ 
let,  Che  primary  fragawncs,  and  Che  stem.  The 
position  (upstream  of  Che  throat)  for  Che  primary 
breakup  moves  from  8.0  mm  to  14.0  am  for  the  ini¬ 
tial  droplet  sizes,  190  pa  to  360  pa,  respectively. 
The  Weber  number  is  baaed  on  the  slip  velocity  and 
initial  droplet  size;  then,  for  the  stem,  Che  Weber 
number  la  baaed  on  an  approximate  frontal  diameter 
(typically  ♦,  *  .30  ♦g).  The  initial  droplets  all 
fail  in  the  Weber  number  range  of  13-20,  and  Che 
sceM  fall  in  a  similar  range.  The  Weber  number 
histories  (Figure  10)  of  all  three  sizes  overlap  to 
within  the  accuracy  of  the  data,  when  the  time  Is 
normalized  by  the  period  of  the  first  natural  fre¬ 
quency  (Eqn.  2).  This  period  for  acoustic  oscilla¬ 
tion  and  Che  breakup  time  are  essentially  equiva¬ 
lent. 

The  fragamnc  size  diitrlbutlon  was  determined 
from  high  magnification  photographs  of  the  recon¬ 
structed  holographic  Images  (ten  events  were 
sampled.).  An  estimate  of  Che  equivalent  spherical 
diameter  %ias  used  for  nonspherical  fragments.  The 
fragoent  number  density  distribution  (Figure  11) 
ranges  from  3-30  pm  and  the  sampling  volume  for 
each  droplet  size  bln  is  corrected  for  Che  depth  of 
field.  The  mean  fragment  size  is  about  10  percent 
of  the  original  droplet  diameter.  Each  fragment 
number  density/  size  distribution  is  a  result  of  a 
single  droplet  breakup  event  in  the  nozzle  contrac¬ 
tion.  The  sample  is  obtained  in  an  instant  in 
time,  but  is  spatially  averaged  along  the  nozzle. 


Droplet  breakup  experiments  have  been  performed 
in  laboratory  scale  nozzle  contraction  (throat  size 
V  1  cm  square).  Room  temperature  air  was  used  for 
Che  test  gas  and  conventional  and  liquid  metal  was 
used  to  produce  the  droplets.  The  droplet  size  was 
varied  from  lUO  pa  to  600  pa  and  the  nozzle  velo¬ 
city  rise  was  varied  from  100  a/s  Co  300  a/s. 
Droplet  velocity  profiles  along  Che  nozzle  were 
resolved  with  laser  velociaetry  and  with  double 
pulsed  laser  holography  which  also  provided  graphic 


images  of  the  droplet  breakup  mechanisms  and  of  the 
droplet  fragnmnts.  The  droplet  velocity  profiles 
were  used  Co  determine  the  Weber  number  profiles 
and  histories.  Thus,  the  critical  Weber  number  was 
determined  and  a  measure  of  the  duration  of  the 
droplet  breakup  event. 

The  breakup  aachaniam  for  conventional  liquid 
droplets  was,  in  general,  equivalent  to  the  Bag 
type  previously  observed  In  shock  Cube  experiments. 
The  critical  Weber  numker  and  breakup  time  data 
revealed  that  when  breakup  occurs  in  a  short  time 
Chat  higher  Weber  numbers  are  required;  and  con¬ 
versely  when  breakup  occurs  in  a  long  time  that 
lower  Weber  numbers  are  required.  This  trend  was 
observed  for  both  water  and  alcohol  droplets.  For 
water  droplets  we  obeerved  a  critical  Weber  number 
of  about  30  for  breakup  times  of  about  2-4  which 
then  decayed  to  a  data  point  for  which  breakup  did 
not  occur  in  the  nozzle  contraction.  One  of  Che 
alcohol  data  points  fell  on  the  water  boundary, 
while  Che  others  reached  Che  highest  critical  Weber 
numbers  observed  in  these  experiMnts.  The  data 
obtained  for  Che  more  viscous  liquid  (a  glycerine/ 
water  mixture)  revealed  the  effect  of  increased 
viscosity  to  be  increased  acceleration.  Viscous 
liquid  droplets  were  observed  to  flatten  to  alisost 
seven  Claes  their  original  dlasmter  before  breakup 
occurred.  The  Increased  acceleration  resulting 

from  Che  larger  cross  section  reduces  the  slip 
velocity  or  Weber  number  achieved  within  the  noz¬ 
zle.  An  additional  effect  of  increased  viscosity 
is  3  reduction  in  the  fragment  size.  The  ratio  of 
the  mean  fragment  to  initial  droplet  size  was 
determined  for  various  Weber  number  failures.  For 
alcohol  and  water  droplets,  Che  ratio  varied  from 
40Z  at  low  Weber  number,  30,  to  about  lOZ  at  high 
Weber  number,  120.  Interesting,  the  high  viscosity 
data  point  reveals  a  fairly  small  fragment  size  for 
Che  low  Weber  number,  40, 

The  liquid  metal  droplet  experiments  required 
higher  dynamic  prsssuras  for  breakup  and  as  a 
result  the  nozzle  was  operated  at  the  sonic  velo¬ 
city.  The  breakup  mechanism  was  observed  to  be  Che 
"Umbrella*  type.  This  eschanlsm  haa  been  observed 
for  mercury  droplets  In  shock  tubs  experiments,  but 
only  for  Weber  numbers  over  80.  The  breakup  mech¬ 
anism  for  lower  Weber  number  shock  Cube  experiments 
"Bag"  type.  A  large  stem  resulted  from  all  primary 
breakup  events  which  eventually  failed  further  down 
Che  nozzle  as  higher  slip  velocities  are  achieved. 
Secondary  breakup  Is  an  leportsnt.  feature  of  the 
SKtal  droplet  breakup  mschanlsm.  In  contrast  to 
Che  conventional  liquid  results,  the  Weber  number 
histories  for  liquid  metals  collapsed  to  a  single 
curve  terminated  by  Che  primary  breakup  event.  All 
three  droplet  sizes  experiaacsd  breakup  at  the  same 
Weber  number,  20.  Also,  cha  duratiom  of  Che  break¬ 
up  event  was  observed  to  ha  a  constant  when  normal¬ 
ized  by  the  natural  period.  The  fragment  size 
distribution  was  observed  to  have  a  mean  value  of 
about  ten  percent  the  original  droplet  size.  Since 
the  critical  Weber  number  did  not  vary  one  might 
expect  the  relative  mean  fragment  size  not  Co  vary. 
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Figure  1.  Nozzle  Facillcy  with  Piezoelectric  Droplet  Injector 


Figure  2.  Holography  Data,  Alcohol  Droplets  (Series  5,  Table  11) 

(a)  Double  Pulse  (it  “  25  us)  Bag  Breakup  Mode 

(b)  Ring  Structure  Observed  Later  in  Breakup  Event 
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Figure  3.  Laser  Doppler  Velocimetry  Data,  Point  Velocity  Distributions  Along 
(See  Figure  1  for  Nozzle  Coordinates) 
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Uater  Droplet  Weber  Number  Profiles 

(Nozzle  Velocity  Rise;  (a)  100  H/S  and  (b)  190  M/S) 
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Figure  5.  Critical  Weber  Numbers  and  Associated  Time  Durations  of  the  Breakup  Event 
Series  Number  (i.e.  1-8)  refers  to  Table  II 


Figure  6.  Mean  Fragment  Size  Dependence  on  Peak  Weber  Number  and  Liquid  Type: 
Open  Symbols  -  LDV  Data;  Closed  Symbols  -  Holography  Data 
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Metal  Droplet  Velocity  and  Weber  Nunber  Profiles; 

(a)  Series  9,  (b)  Series  10.  and  (c)  Series  11  (Table  II) 
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